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Abstract 

Two recently proposed expressions for the computation of the entropy in the microcanonical 
ensemble are compared, and their equivalence is proved. These expressions are valid for a certain 
class of statistical mechanics systems, that can be called mean-field-like systems. Among these, 
this work considers only the systems with the most usual hamiltonian structure, given by a kinetic 
energy term plus interaction terms depending only on the configurational coordinates. 
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The micro canonical ensemble generally is not the chosen ensemble for the computation 
of thermodynamic observables of statistical systems. Exact or approximate analytical pro- 
cedures are usually easier in the canonical ensemble. However, a recent surge of interest 
for the statistical properties of systems with "unconventional" thermodynamics has made 
desirable to perform calculations in both ensembles. These systems are characterized by the 
fact that they are small ll, and thus surface effects cannot be neglected, or by the long 
range nature of the interaction between their microscopic cases 
additivity of the energy and of other thermodynamic functions does not hold. 

The most relevant feature of long range systems is that very often, even in the thermody- 
namic limit, canonical and microcanonical ensembles are not equivalent, and this justifies the 
necessity to perform calculations in both ensembles: contrary to the usual case, the system 
can have different values of the macroscopic observables, if it is considered at fixed energy 
or at fixed temperature. The issue of ensemble inequivalence, first discussed theoretically by 
Hertel and Thirring ^ | , has later been thoroughly treated in the mathematical physics lit- 
erature (see Refs. y] and references therein). The inequivalence has been put in evidence 



both in self-gravitating systems 



nij, 



y, Q, ^ and in other, magnetic- type, systems j^i 0, 
which are often easier to study. A fundamental tool for the mathematical analysis of en- 
semble inequivalence has been found in the theory of large deviations ^Uj, as applied to 
statistical mechanics (see Ref. and references therein): this application employs the 
important concept of macrostates of the systems 0, 14, 1^ . 

In simple terms, nonequivalence can be ascribed to the nonconcavities that are found in 
the microcanonical entropy |4]. Not only does this lead to different results for the entropy, 
if computed in the microcanonical and in the canonical ensembles, but it can also be the 



cause, in the microcanonical ensemble, of negative specific heats l6|,l7|,|8|,|9|,lll|,ll6|,|l3,Q. It 

has also been found that ensemble inequivalence is associated to the presence of first order 
phase transitions in the canonical ensemble 



Through the use of large deviation techniques, it is possible to derive a general expression 
for the microcanonical entropy of statistical mechanics systems belonging to a certain class [J, 



IG 



[fsj ]. that we can call mean-field-like systems. This terminology is motivated by the 
observation that in these systems the interactions can be thought of as giving rise to a 
mean field acting on the particles, although this may not be easily visible as in, e.g., some 
simple lattice magnetic models. If one accepts to sacrifice a full mathematical rigour, the 
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derivation of this expression can be short jlO|, ll9| • Recently, another expression has been 
proposed j^], starting from the form that the canonical partition function assumes in mean- 
field-like systems; no obvious and self-evident relation with the large deviation expression 
can be easily found. In this work our purpose is to make a connection between these two 
different expressions. Our motivation is to show the equivalence of the two approaches, and 
at the same time to give an operative procedure for the derivation of the particular form 
of the canonical partition function in mean-field-like systems, that in Ref. j20| is used as a 
starting point, but is not derived. 

In a system with degrees of freedom and canonical coordinates {{qi,Pi)}, we are inter- 
ested in the computation of the entropy per particle in the thermodynamic limit: 

.(e) = lim 1 In / {Ul,dq,dp,) 6 {H{{{q,,p,)}) - Ne) , (1) 

where H is the Hamiltonian and e is the energy per particle (or energy density). Let us first 
define the class of systems that are considered here. We suppose that the Hamiltonian of 
our systems can be written in the form: 

n 

H = NY,9k{Mk), (2) 

k=l 

with 

1 ^ 

Mfc = -5^/fc(g„p,), (3) 

i=l 

where fk and Qk are smooth functions. It should be pointed out that the theory of large 
deviations can be applied to systems defined by looser conditions on the Hamiltonian; in 
particular, the form Q could also be the leading term of an expansion in A^, in which 
the relative weight of the remaining terms vanishes in the thermodynamic limit; even the 

n 

simple additive structure displayed in (j2I) is not necessary 4]. However, for simplicity here 
we restrict to systems in which Eq. Q holds exactly; moreover, we treat the case in which 
the k = 1 term is the kinetic energy, i.e., fi{q,p) = flip) = (with a suitable choice 
of the mass unit) and gi{x) = x, while the functions fk for A; > 2 represent the interaction 
terms, and depend only of the configurational coordinate q. Thus, we are considering mean- 
field-like systems in which each particle is described by one degree of freedom. Later it will 
be shown that the extension to more than a degree of freedom per particle (e.g., cartesian 
coordinates of particles in ordinary threedimensional space) is easily performed. 
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We now show our first expression, that we call the large deviation expression of the 
entropy per particle. The cited literature can be consulted for its derivation. One finds 
that: 



where 



and where 



s(e) = sup J(mi, . . . ,m„), (4) 

mi,...,m.n,T, 9k{-mk)=e 



J{mi,...,Tnn) = inf L(zi, . . . , z^, mi, . . . , m„) (5) 



n 



. . . ,2;„,mi, . . . ,m„) = ln(e^fe=i^'=-^''-)i - \^^^Zkmkj; (6) 
here is the one particle phase space (unnormalized) average defined by: 

{■),= fdqdpi-). (7) 



Of course, depending on the form of the functions fk, there can be restrictions on the sign 
of the ZfcS, in order to have convergence of the one particle phase space average in the right 
hand side of Eq. (jH)) . 

The second expression is obtained, as already mentioned, in Ref. [2^, and it is derived 
under the hypothesis that the canonical partition function of the system can be written in 
the form: 

Z{P) ^ I {Ul,dq,dpi) exp[-PHi{iq,,pi)})] (8) 

= C j dyi...dyrexp[N(j){p,yi,...,yr)], 

where the factor C has at most a power law dependence on A^, so that it does not contribute, 
in the thermodynamic limit, to the free energy per particle f{f3), given by: 

-(3f{(3) = hm ^lnZ{(3) 

= sup (l){(3,yi,...,yr). (9) 

In the derivation of the microcanonical entropy, apart from the assumption of analiticity, 
there is no other hypothesis about the form of the real function 0, about the number r of 
"order parameters", and about the way in which has been obtained. Actually, Ref. 



treats the case r = 1, but the generalization to r > 1 is straightforward 



2l\ . It can then be 
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shown that the microcanonical entropy per particle is given, in the thermodynamic hmit, 
by: 



sup 
yi,---,yr 



inf {/3e + (f){/3,yi,...,yr)) 



(10) 



In Ref. the systems for which Eq. (jSJ holds are called mean-field-like systems; thus we 
will call Eq. (fTUI) the mean-field-like expression of the entropy per particle. Our purpose 
here is to show that this last expression is equivalent to Eqs. (jHIHl)- 

First of all, let us specialize the expressions to our case, that is concerned with the 
most usual situation in which in the Hamiltonian there is a kinetic energy term plus terms 
depending only on the configurational coordinates and giving the interaction. Then, in the 
one particle phase space average of Eq. the integrals in q and p separate; besides, the 
minimization problem defined in Eq. (0) decouples in one concerning zi, that is easily solved, 
and one concerning Z2, ■ ■ ■ ,Zn- One ends up with: 



J(mi, . . . ,m„) 



+ inf 

Z2,...,Zn 



= ^ln(27re) + ^ln(2mi) 
G{Z2, ...,Zn,m2,.. .,mn), 



(11) 



where 



G{z2, . . . ,Zn,m2, . . . , rUn) = G{{zk}, {rrik}) 



Zknik 



(12) 



,k=2 



where now is the analogous of Eq. ((7j), but with the integration only in q. Let us denote 
with 'Zk, k = 2, . . . ,n, the solution of the minimization problem defined in the right hand 
side of Eq. (jlip : the ^^s are thus functions of m2, . . . yfrin- Furthermore, since gi{x) = x in 
Eq. (j2I), the maximation problem in Eq. can be performed by a free maximation over 
m2, . . . , TTin after putting: 

n 

mi = e - ^5(fc(mfc). (13) 

k=2 

Therefore, neglecting the constant term in the right hand side of Eq. (fTTjl . we obtain: 



s(e) = sup 

m2,...,m„ 



1 " 

-In (2e - 2^5(fc(mfc) 

k=2 



+ G{z2,...,Zn,m2, 



(14) 



5 



This is the large deviation expression of the thermodynamic hmit of the entropy per particle, 
for the class of systems considered here. 

We can now show, starting from the definition of the canonical partition function, that 
Eq. (jlUj) can be transformed in Eq. (jl4|) : moreover, we can derive a definite form of (f) in 
Eq. (jS)). In the computation of the canonical partition function, it is possible to adopt a 
procedure that is similar to that employed in the derivation of the large deviation expression 
of the entropy In fact, taking into account the form of the Hamiltonian, Eqs. (j2I) and 
dni), and specializing to our case, we can write: 

Z{[3) = f ^ ] [ dgi . . . dgArdm2 . . . dm„ (15) 

n 



{lll^,N6{NMk{{q,}) - Nmk)) exp 



k=2 



The phase space coordinates are now only in the 6 functions. Using the inverse Laplace 
transform for these functions we find that: 



dgi . . . dqN{Ul^2N5iNM,i{q,}) - iVm^)) 

= J^dz2... dzn (^^^ exp [NGiizk}, {m,})] , (16) 

with the function G defined in fll2j) . and where the integration path for each of the variables 
Zk is the imaginary axis, or a path obtained by deformation of the imaginary axis without 
crossing singularities. We compute the integral in the right hand side of for —>■ oo, 
with the saddle point method (full mathematical rigour would require to prove the analiticity 
of the function G). To avoid unphysical oscillatory behavior, the relevant stationary point 
must lie on the real Zk axes. It can then be inferred that this stationary point can be found 
looking for the minimum of G with respect to real values of Z2, . . . , Zn- This can be deduced 
in the following way. Let us first define the weighted one particle average (■)^ by: 

(.) = Alii 111 (17) 

\/UI I ^ f \ • \ ' J 

The leading order variation AG around a stationary point of the function G is easily found, 
from Eq. (jl2|) . to be: 



n n 



AG = }^}_^{{fk - {fkU ifi - {fiU)n.dzkdzi 

k=2 1=2 



Restricting to real ZkS we have both that G is a real function, and that Eq. ()18|) is a positive 
quantity, meaning that a stationary point is a minimum. Thus the saddle point evaluation of 
the integral in the right hand side of Eq. is obtained traversing the relevant stationary 
point along a line parallel to the imaginary axis for each z^, so that the real part of G 
attains the maximum in the integration domain. At the end, keeping only the exponential 
dependence on N, Eq. (fT^ becomes: 

y dgi . . . dqN{lll^2mNMk{{q,}) - Nnik)) 
oc exp [NG(z2, . . . ,Zn,m2, . . . , m„)] . (19) 

This result can be inserted in Eq. (I15j) . and one therefore obtains an expression like in 
Eq. (jS)), with r = n — 1 and, after the change of notation yi, . . . ,yr 1712, ... , m„, with 
given by: 

0(/?,m2,...,m„) = -In ( — j -/3^^(mfc) 

^ ^ ^ k=2 

+ G{Z2, . . . ,Zn,m2, . . . ,mn). (20) 

At this point, with the very simple dependence on /3 of this expression, it is straightforward 
to perform the minimization problem in (3 indicated in Eq. (fTUI) : we obtain: 

1 " 

= 2e-2Y,9imk). (21) 

^ k=2 

This is substituted in (3e + (f){(3, m2, . . . , rrin), and then the mean-field- like expression of the 
microcanonical entropy per particle, Eq. (jlOp . becomes, again neglecting the same constant 
term as before, exactly the same as the large deviation expression, Eq. (fT^ . 

The extension to systems in which the particles move in a D-dimensional space is straight- 
forward. In this case the canonical coordinates {qi,Pi} in Eqs. (HJ), (jS)), and (jHI) are to be 
intended as D-dimensional vectors {qj,Pi}. Analogous replacements have to be done for 
the configurational coordinates qi in Eq. p6|) . in the one particle phase space average 
(see Eq. ((7j)), and in its configurational restriction {■)iq. With e we still denote the energy 
per particle (not per degree of freedom). At the end, the common result for the large de- 
viation and the mean-field-like expressions of the entropy per particle, Eq. p4|l . is simply 
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generalized, neglecting constant terms, to: 



s(e) = sup 

m2,...,m„ 



D 
2" 



In (2e - 2^gk{mk) 

k=2 

+ G{z2,...,Zn,m2,...,mn) 



(22) 



We end with the following observation. The form obtained for the function in Eq. ()20|) 
is useful for our proof of equivalence of the two approaches to the computation of the mi- 
crocanonical entropy per particle, but by no means it is unique. For example, with simple 
quadratic functions gk in Eq. 0, it is common to employ the Hubbard- Stratonovich trans- 
formation [2^ to linearize the dependence of the Hamiltonian on the phase space variables. 
One then arrives, for the canonical partition function, to an expression as in the righmost 
hand side of Eq. (jSJ, with the functional form of 0, and the physical meaning of the "order 
parameters" mfc, that can be different from that in Eq. (jlOJ. Depending on the concrete 
cases, the function cj) obtained may result in an expression that is easier to analyze numeri- 
cally than Eq. 
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